Flow past superhydrophobic surfaces with cosine variation in local slip length 
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Anisotropic super-hydrophobic surfaces have the potential to greatly reduce drag and enhance 
mixing phenomena in microfluidic devices. Recent work has focused mostly on cases of super- 
hydrophobic stripes. Here, we analyze a relevant situation of cosine variation of the local shp length. 
We derive approximate formulae for maximal (longitudinal) and minimal (transverse) directional 
effective slip lengths that are in good agreement with the exact numerical solution and lattice- 
Bolzmann simulations for any surface slip fraction. The cosine texture can provide a very large 
effective (forward) slip, but it was found to be less efficient in generating a transverse flow as 
compared to super-hydrophobic stripes. 



I. INTRODUCTION 

The design and fabrication of textured superhydropho- 
bic surfaces have received much attention in recent years. 
If the recessed regions of the texture are filled with gas 
(the Cassie state) , roughness can produce remarkable hq- 
uid mobihty, dramaticaUy lowering the abihty of drops 
to stick [T]. These surfaces are known to be self-cleaning 
and show low adhesive forces. In addition to the self- 
cleaning effect, they also exhibit drag reduction for fluid 
flow. Thus, they are of importantance in the context of 
transport phenomena and fluid dynamics as well [21 [3 H] . 
Many sea animals e.g. shark and other flsh are known to 
possess superhydrophobic skin 5j. Also many artificial 
textures have been designed to increase drag reduction 
efficiency [6]. 

This drag reduction is associated with the liquid slip- 
page past solid surfaces. This slippage occurs at smooth 
hydrophobic surfaces and can be described by the bound- 
ary condition [21 13 IH|, "slip = bdu/dz, where Usu-p is 
the (tangential) slip velocity at the wall, du/dz the lo- 
cal shear rate, and b the slip length. A mechanism for 
hydrophobic slippage involves a lubricating gas layer of 
thickness 5 with viscosity /ig much smaller than that of 
the liquid ^i 0, so that b ~ 5{^l/^JLg - 1) ~ 505 HH]- 
However, at smooth flat hydrophobic surfaces 5 is small, 
so that b cannot exceed a few tens of nm [TTl |T2 [13 E] . 
In case of superhydrophobic surfaces, the situation can 
change dramatically, and slip lengths up to tens of /im 
may be obtained over a thick gas layer stabilized with a 
rough texture [HI [16] . 

To quantify the flow past heterogeneous surfaces it 
is convenient to apply the concept of an effective slip 
boundary condition at the imaginary smooth homoge- 
neous, but generally anisotropic surface [21 [II]. Such an 



effective condition mimics the actual one along the true 
heterogeneous surface, and fully characterizes the real 
flow. The quantitative understanding of the effective slip 
length of the superhydrophobic surface, boff , is still chal- 
lenging since the composite nature of the texture in ad- 
dition to liquid-gas areas requires regions of lower local 
slip (or no slip) in direct contact with the liquid. For an 
anisotropic texture, the effective slip generally depends 
on the direction of the flow and is a tensor, bos = {bf^} 
represented by a symmetric, positive definite 2x2 ma- 
trix [18] 

diagonalized by a rotation 




Eq.([T]) allows us to calculate an effective slip in any di- 
rection given by an angle 9, provided the two eigenvalues 
of the slip- length tensor, b\^ (6 ~ 0) and 6^ [Q = 7r/2), 
are known. The concept of an effective slip length tensor 
is general and can be applied for an arbitrary channel 
thickness [IH], being a global characteristic of a chan- 
nel [3] , so that the eigenvalues normally depend not only 
on the parameters of the heterogeneous surfaces, but also 
on the channel thickness. However, for a thick (compared 
to a texture period, L) channel we are interested in here 
they become a characteristics of a heterogeneous inter- 
face solely. 

In case of an anisotropic isolated surface (or a thick 
channel limit) with a scalar local slip b{y), varying in 
only one direction, the transverse component of the slip- 
length tensor was proven to be equal to a half of the 



2 







h" bi 






^ — j#///IlllV ^ 




















y 

> 



FIG. 1: (Color online) Sketch of the SH surface with a co- 
sine relief, and its equivalent representation in terms of flow 
boundary conditions. 



longitudinal one with twice larger local slip, 2b{y) [20] 



(2) 



A remarkable corollary of this relation is that the flow 
along any direction of the one-dimensional surface can 
be easily determined, once the longitudinal component of 
the effective slip tensor is found from the known spatially 
nonuniform scalar slip. 

One-dimensional superhydrophobic surfaces are very 
important for a class of phenomena, which involves 
"transverse" hydrodynamic couplings, where an applied 
pressure gradient or shear rate in one direction gener- 
ates flow in a different direction, with a nonzero per- 
pendicular component. This can be used to mix adja- 
cent streams, control the dispersion of plugs, and posi- 
tion streams within the cross section of the channel [21] . 
Such grooved surfaces can be easily prepared by modern 
lithographic methods 6]. 

Most of the prior work focussed on a flat, periodic, 
striped superhydrophobic surface, which corresponds to 
patterns of rectangular grooves. The flow past such 
stripes was tackled theoretically [3t [22 l l23 l l2t | [25 l l26] . 
and several numerical approaches have also been used 
either at the molecular scale, using molecular dynam- 
ics ^27], or at larger mesoscopic scales using finite ele- 
ment methods [351 [2S] , lattice Boltzmann J_9 and dissi- 
pative particle dynamics [30] simulations. For a pattern 
composed of no-slip (5 — 0) and perfect-slip {b — oo) 
stripes, the expression for the eigenvalues of the effective 
slip-length tensor takes its maximum possible value and 
reads [H [ST] 



^dcal — 2^idoal - 



sec 



7r(l 



(3) 



and fo[ijeai to depend only logarithmically on (j> and scale 
as —Lint/). At a qualitative level, this result means that 
the effective slip lengths essentially saturate at the value 
fixed by the period of the roughness. In case of stripes the 
perturbation of piecewise constant local slip has a step- 
like jump on the heterogeneity boundary, which leads to a 
singularity both in pressure and velocity gradient |20| by 
introducing an additional mechanism for a dissipation. It 
is natural to assume that an anisotropic one-dimensional 
texture with a continuous local slip could potentially lead 
to a larger effective tensorial slip. 

In this paper we address the issue of the effective 
slip of flat surfaces with cosine variation in the local 
slip length, which corresponds to modulated hydropho- 
bic grooved surfaces with a trapped gas layer (the Cassie 
state) as shown in Fig. nl Flows over hydrophilic sur- 
faces (the Wenzel state) with cosine surface relief of 
small amplitude have been studied by a number of au- 
thors [m Ha [33 IMl ESI- Previous studies of similar 
grooves in the Cassie state have investigated only small 
variations in local slip length [351 [35] . We are unaware of 
any previous work that has studied the most interesting 
case of flnite and large variations in the amplitude of a 
local cosine slip. 



II. THEORY 

Consider a shear flow over a textured flat slipping 
plate, characterized by a slip length b{y), spatially vary- 
ing in one direction, and the texture varying over a period 
L (as shown in Fig. T]). We use a rectangular coordinate 
system (x, y, z) with origin at the wall. The z— axis 
is perpendicular to the plate. Our analysis is based on 
the limit of a thick channel or a single interface, so that 
the velocity profile sufficiently far above the surface, at 
a height large compared to L, may be considered as a 
linear shear flow. All variables are non-dimensionalized 
using the texture period L as the characteristic length, 
the shear rate of the undisturbed flow G and the fluid 
viscosity ^. 

The dimensionless fluid velocity is sought in the form 

V = U + Usiip + ui {x, y, z) , 

where U — zei, I = x,y is the undisturbed linear shear 
flow, and e; are the unit vectors parallel to the plate. The 
perturbation of the flow, which is caused by the presence 
of the texture, involves a constant slip velocity Ugiip = 
(wsiip, Wgiip, 0) and a varying part Ui — {u,v,w) of the 
velocity field. A periodic velocity Ui should decay at 
infinity and has zero average: 



Uidy = 0. 



(4) 



where (j) is the fraction of the no-slip interface. In the 
limit of vanishing solid fraction, it therefore predicts b^[^^^\ 



At a small Reynolds number Re — GL^/i', Ui satisfies 
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the dimensionless Stokes equations, 



A. Longitudinal configuration 



V • ui - 0, 
Vp - Aui = 0. 



(5) 



The boundary conditions at the wall and at infinity are 
defined in the usual way as 



z = G: Usiip + \iir~(i (y) 
w = 0, 



duir 

dz 



ui = 0, 



P{y)eu (6) 
(7) 



(8) 



where Uit- = {u, v, 0) is the velocity along the wall and 
(3 = b/L is the normalized local slip length. The eigen- 
values of the effective slip-length tensor can be obtained 
as the components of Ugiip : 



Lv, 



slip ■ 



(9) 



III. COSINE SLIP LENGTH 



In this Section we consider a ID periodic texture with 
the local slip length 



/3 = /3o + 2/3i cos (27ry) . 



(10) 



The coefficients should satisfy /?o > 2/3i > 0, in order to 
obey /3 (y) > for any y. 

The disturbance velocity field is presented in terms of 
Fourier series as 

oo 

Ui = u* (n, z) exp (z27rny) , u* = {u* ,v* ,w*) . 

n— — oo.n^O 

(11) 

A general solution of the Stokes equations for the longi- 
tudinal flow, U = zBx decaying at infinity, reads ^20] 



u* = Xn exp (— 27r \n\ z) , 



(12) 



V* =w* ^ 0, 

and that for the transverse flow, U = zSy is given by 

u* = 0, 

u* = r„ exp (-27r |n| z) (1 - 27mz) , (13) 

w* = — i27rnF„z exp (— 27r |n| z) . (14) 

The Fourier coefficients X„ and Yn are determined from 
the Navier slip boundary condition (|6|. 



Since the local slip length is an even function of y, 
the solution (111 is also an even function. This requires 
Xn = X-n, SO it is sufficient to evaluate X„ for n > 0. 
The Navier slip boundary condition ^ can be written, 
following to [5D] , in terms of Fourier coefficients as a lin- 
ear system 



"slip = ^0 - 2eiXi, 



diXi + 62X2 = /3i 



(15) 



(16) 



n > 1 : e„_iX„_i + d„X„ + e„+iX„+i = 0, (17) 



1 + 27rn^o, 



27rn^i. 



(18) 



A three-diagonal infinite linear system ( 15 1-( 17 1 should 



be solved numerically to find teh unknown A"„ by trun- 
cating the system. 

In the limit of large slip, /3q > 2/3i ^ 1, the asymp- 
totic solution to (15l-(17) can also be constructed. To 



the leading order in , the first term in ( 18 1 can be 
neglected compared to the second one, and the system 
(16l-(17) is rewritten for new variables t„ = 2TmXn as 



h + Xt2 = A, 



n> 1 : Xtn-l +tn+ Xtn+l 



0, 



(19) 



(20) 



where A = /3i//3o < 1/2. The solution of the last system 
is a geometric progression, t„^i = q"ti, with 



Q = 



-A-i + VA-2 -4 



A 



1 + qX Po + qpi 

Therefore, the final expression for the slip length, = 
iwsiip, in view of (15), takes the following form at /3o > 
2/3i > 1 : 



4bl 



(21) 



The asymptotic solution for the velocity field can also 



be derived from Eqs.(ll| and (12). The normal compo- 



nent of velocity gradient can be written as 



du 
dz 



= -2ti q"-^Re {exp [27rn {iy - z)]} 



(22) 



n=l 



The factor 2 is due to the contribution of n < 0. The sum 



in ( 22 ) is also a geometric progression, so that 

du 2ti exp (— 27rz) [cos (27ry) — q exp (— 27rz) 
dz s 



, (23) 
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s = 1 — 2g cos {2TTy) exp {—2nz) + exp {—At:z) 
The last equation can be integrated over z to give 



27rg 



■ In ; 



(24) 



The asymptotic solutions (21), (23) and (24) predict well 



the numerical data even at finite 6o (see next Section). 



B. Transverse configuration 

It was found by [20! that the velocity components for 
the transverse flow can be expressed in terms of the 
longitudinal one calculated for twice larger local slip, 
U2 = u [2/3 (y)] : 



1 / du^ 
2^ + ^97 
z du2 



Using ^ and ( 25 ) we derive 



(25) 

(26) 
(27) 

(28) 



Thus, the texture is isotropic at large h^. One can demon- 
strate that the conclusion about isotropy of the slip- 
length tensor is general and valid for any textures (see 
Appendix [A| . 

The values A, g, <i remain the same for the double slip 
length since they depend on the ratio /3i//3o only. As a 
result, we obtain 

U2 — u — — — In s, 
2tt 



in s 

zti exp (— 27rz) [cos (27ry) — (/exp (— 27rz)] 
s 

zti exp (— 27rz) sin {^ny) 



w 



IV. SIMULATION METHOD 



(29) 



(30) 



For the modelling of fluid flow in a system of two paral- 
lel plates, we employ the lattice Boltzmann (LB) method 
[37] . Lattice Boltzmann methods are derived by a phase 
space discretization of the kinetic Boltzmann equation 



' d_ 
di 



V- Vr 



/(r,v,i)-f^, 



(31) 



which expresses the dynamics of the single particle prob- 
ability density f{r,v,t). Therein, r is the position, v 
the velocity, and t the time. The left-hand side models 
the propagation of particles in phase space, the collision 
operator on the right hand side accounts for particle 
interactions. 

Constructing the lattice Boltzmann equation, the time 
t, the position r, and the velocity v are discretized. This 
discrete variant of Eq. (31 ) 



/fe(r + Cfe, t + 1) - A(r, t) - rife. A: - 0, 1, ... , B, 

(32) 

describes the kinetics in discrete time- (At) and space- 
units (Ax). We employ a widely used three-dimensional 
lattice with B = 18 discrete non-zero velocities (D3Q19) 
which is chosen to carry sufficient symmetry to allow for a 
second order accurate solution of the Navier-Stokes equa- 
tions. Here, for ft, we choose the Bhatnagar-Gross-Krook 
(BGK) collision operator [38] 



T 



(33) 



which assumes relaxation on a linear timescale r towards 
a discretized local Maxwell-Boltzmann distribution /^'. 



The kinematic viscosity 



of the fluid is related to 



the relaxation time scale. In this study it is kept constant 
at r = 1.0. 

Stochastic moments of / can be related to physical 
properties of the modelled fluid. Here, conserved quan- 
tities, like the fluid density p{r,t) — Po X^fe /*:(i"j ^^^d 
momentum p(r, i)u(r, i) = Po Sfc Cfe/fe(r, i), with po be- 
ing a reference density, are of special interest. 

Slip over hydrophobic surfaces is commonly modelled 
by introduction of a phenomenological repulsive force 
[Ml EOl mi m ig. The magnitude of interactions 
between different components and surfaces, as deter- 
mined by simulation parameters allows to specify arbi- 
trary contact-angles [31113111151135]. Other approaches 
include boundary conditions taking into account specular 
reflections A7, 48, 49J, or diffuse scattering jSO^^l, 52j, 
respectively. The strategy applied for this work employs 
a second order accurate on-site fixed velocity boundary 
condition to simulate wall slippage. Here, the velocity at 
the boundary is set proportional to the local stress im- 
posed by the flow field as well as a slip length parameter 
adjusting the stress-response. For the details of the im- 
plementation we refer the reader to [SJl [S3] • Local slip 



lengths are calculated according to Eq. ( 10 ). Varying slip 
patterns are applied to the x-y-plane a,t z = 0. Periodic 
boundary conditions are employed in x and y-direction, 
thus reducing the simulation domain to a pseudo-2D sys- 
tem. By exploiting the periodic boundaries, only one 
single period needs to be resolved. The Couette fiow is 
driven by applying a constant velocity of v = 0.1 (in 
lattice units) in the x-y-plane at z = Zmax- 

The resolution of the simulated system is given by the 
lattice constant 



Ax 



H 



(34) 



5 




1.5r 



0.1 0.2 0.3 0.4 

b /b 

10 



FIG. 2: Eigenvalues of the effective slip-lengtfi tensor as a 
function of bi simulated at fixed 6o = 1 (symbols). The lon- 
gitudinal effective slip length, fe|!g, is shown by circles, and 
the transverse effective slip, b^g is presented by diamonds. 
Solid and dashed curves denote the corresponding theoretical 
values obtained by numerical Fourier-series solutions. The 
asymptotic (isotropic) solution, Eq. (281, expected in the 



limit bo 3> i is shown by the dash-dotted line. 



where is the number of discretization points used 
to resolve the height of the channel. As we consider 
thick channels, we choose a height to periodicity ratio 
of H/L = 10. 

The number of timesteps required to reach a steady 
state depends on the channel height, the velocity of the 
flow as determined by the driving acceleration as well as 
the fraction of slip and no slip area at the surface. We 
find that in order to reduce the deviation from theoretical 
predicted values below 5 percent, a domain size of 96 x 
1 X 960Aa;^ is required. For this geometrical setup, with 
a shear rate in the order of 7 = 1 • 10^** (in lattice units), 
a simulation time in the order of 10 million timesteps 
is needed to reach the steady state. We compare the 
theoretical predictions by measurements of 6eff obtained 
from velocity profiles. The profile of the linear shear 
velocity is fit by a linear function in the region far from 
the surface pattern. From this fit the effective slip lengths 
and shear rate are determined. 



V. RESULTS AND DISCUSSION 

In this section, we present the LB simulation results 
and compare them with predictions of the continuous 
theory. 

We start with varying the amplitude of cosine pertur- 
bations of the shp length, 61, at fixed 60 = 1. Fig. [2] 
shows simulation data for b^g and b^g as a function of 
bi/bo. These results show that the largest possible value 
of bcs/bo is attained when 61 = 0, i.e. for a smooth 
hydrophobic surface with b{y) = bg. In this situation 
the effective slip is (obviously) isotropic and equal to the 
area-averaged slip &o- When increasing the amplitude 61, 
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FIG. 3: (a) Effective slip lengths computed for 61/60 = 0.5, 
which correspond to a texture with no-slip lines, vs. average 
slip. The notations are the same as in Fig. [2] Dash-dotted 
and dotted lines show the effective lengths for longitudinal 
and transverse stripes as a function of 1/(30) calculated with 
Eq.Q. (6) The cosine profile of the local slip length with 
60/-L — 5, 61/60 — 0.5 (solid curve) and the stripe profile with 
(j) = 0.06 (dashed line) with the same longitudinal effective 
slip lengths. 



there is a small anisotropy of the flow, and the eigenvalues 
of the slip-length tensor decrease. Therefore, in the pres- 
ence of a cosine variation in slip length the effective slip 
always becomes smaller than average. This conclusion 
is consistent with earlier observations made for different 
textures [31 [5S]. To obtain theoretical values, the lin- 



ear system, Eqs. (15)-(17l has been solved numerically 
by using the IMSL-DLSLTR routine. We sec that the 
agreement is excellent for all 61/feo, indicating that our 
asymptotic theory is extremely accurate, and confirming 
the relation ( 25 1 between the longitudinal and transverse 
slip lengths. Also included in Fig.[2]is the asymptotic for- 
mula ( 28 ) obtained in the limit of large bo. Note that this 



formula is surprisingly accurate even in the case of finite 
60, except for the texture with &1/60 — 1/2 (no-slip point 
at y = 1/2). 

Fig. [3] (a) shows the simulation data for effective slip 
lengths as a function of average slip, 60/^1 for a tex- 
ture with the no-slip point {bi/bo — 1/2). Also included 
are theoretical (Fourier series) curves. The fits are quite 
good for bo/L up to 10, but at larger average slip there 



6 



is some discrepancy. The simulation results for 6"g and 
b^ff give smaller values than predicted by the theory. A 
possible explanation for this discrepancy is that the ma- 
jor contribution to the shear stress at large bo/L and 
bi/bo = 1/2 comes from a very small region near the no- 
slip point (as we discuss below). The discretization error 
of the LB simulation becomes maximal in this region, 
and is particularly pronounced for the velocity gradients 
of systems with large effective slip. While we observe 
deviations around the no slip extremal value, the curves 
converge fast when stepping away from it and the ex- 
cellent agreement of the measured effective slip suggest 
that the influenece of discretization errors on the mean 
flow is negligible at the resolution used. The asymp- 
totic formula, Eq.(28), predicts bes = 0. This likely in- 




dicates that in this situation it is necessary to construct 
the second-order term of expansions for eigenvalues. At 
relatively large bo/L, the effective slip lengths can be well 
fitted as 

blg/L ~ 0.1871 + 0.3175 In (6o/i + 1.166) , (35) 
6fff/L ~ 0. 2036 -I- 0.158 8 In (6o/L -1-0. 583) . 

In other words, they scale as \n{bo/L) at large bo/L. 

Since the effective slip lengths for a texture decorated 
with perfect-slip stripes, Eq.([3]), also shows a logarith- 
mic growth (with (j)), in order to compare these two one- 
dimensional anisotropic textures the theoretical curve for 
stripes is included in Fig.jsja). It can be seen that in the 
limit of large average slip the asymptotic curves for longi- 
tudinal effective slip for stripes and cosine texture nearly 
coincide. This means that both textures generate the 
same forward flow in the longitudinal direction. Simple 

estimates suggest blf^ (/3o) ^ ^I'doai [1/ i^M] ■ Perhaps the and ^ are in a good agreement with numerical results 
most interesting and important aspect of this observation 
is that, from the point of view of the longitudinal effec- 
tive slip, the "wide" cosine texture with f3o — 5 taken for 
our numerical example is equivalent to the patterns of 
stripes with the extremely low fraction of no-slip regions, 
</> = 0.06 (see Fig.|3];5)). These results may guide the de- 
sign of superhydrophobic surfaces for large forward flows 
in microfluidic devices. Note, however, that in the situa- 
tion when longitudinal slip for both textures are similar, 
the cosine texture shows a larger transverse effective slip 
as seen in Fig. ^a) . This means that textures with the 
cosine variation in the local slip length are less anisotropic 



FIG. 4: (a) The velocities and (6) the normal velocity 
gradients along the wall for the textures with bi/bg = 
1/3, bo/L = 0.2 (solid curve, diamonds), bo/L — 1 (dashed 
curve, circles), bo/L = 5 (dash-dotted line, crosses). Do tted 
curves show predictions of asymptotic formulae, Eqs. ( 24 1 and 
(23|. 



than stripes, and blg/b^g < &Lai/feMoai = 2- Therefore, 
cosine textures are less optimal for a generation of robust 
transverse flows as compared with a sharp-edge stripe ge- 
ometry. 

The flow direction is associated with hydrodynamic 
pressures in the film, which is related to the heteroge- 
neous slippage at the wall. Fig. [i] shows the profiles of 
the velocity and of the normal velocity gradient along 
the wall for different (3o and /3i//3o > 1/2- The veloc- 
ity dependence u{x,y,0) is smooth, and j^{x,y,0) is 
finite for any /3q and unlike the striped textures with 
piecewise-constant /3 [5D]. Asymptotic predictions (24) 



and simulation data. 

Similar theoretical and simulation results, but ob- 
tained for a texture with no-slip point, Pi/Po = 1/2, 
are shown in Fig. [5] In this situation we find that 
du/dz (1/2) = 27r^o for all Pq. 

Finally, we should like to stress that a very small region 
near the no-slip point gives a main contribution to the 
shear stress at large b^/L. For the major portion of the 
texture far from this region, we have du/dz ~ —1, so 
that the total shear stress is zero, and this part of the 
texture is shear-free. Since the maximum values of the 
normal velocity gradients grow like b^/L one can expect 
that a length scale of this small region is L^/bo ^ L, or, 
equivalently, the curvature radius, r — {cPb/dy^^ — 
L"^/ (47r^&o) ^ L, at the no-slip point. The validity of 
this assumption is justified in Fig. [6j where the gradients 
for several values of bo/L versus stretched coordinates, 
bo {y ~ 0.5) /L, are presented. The curves are very close 
for bo/L > 5. Therefore, the normal gradient distribution 
and the dimensionless effective slip length in this case are 
controlled by the ratio r/L only. These conclusions can 
be extended to any b {y) characterized by a small radius 
r = (^cPb/dy^^ <^ L near the no-slip point and by a 
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y 

FIG. 5: (a) The velocities and (fo) the normal velocity gra- 
dients along the wall for the no-slip textures (61/60 = 1/2). 
Other notations are the same as in Fig. [4] 



slip length, and have evaluated resulting effective slip- 
page and the flow velocity. We have found that the 
cosine texture can provide a very large effective (for- 
ward) slip, but generates a smaller transverse velocity 
to the main (forward) flow than discrete stripes consid- 
ered earlier. Our approximate formulae for longitudinal 
and transversal directional effective slip lengths are val- 
idated by means of lattice Boltzmann simulations. Ex- 
cellent quantitative agreement is found for the effective 
slippage as well as for the flow fleld. Slight deviations 
of the observed velocity gradient close to the no-slip ex- 
tremal value can be explained by discretization errors. 

Acknowledgments 

This research was partly supported by the Russian 
Academy of Science (RAS) through its priority program 
'Assembly and Investigation of Macromolecular Struc- 
tures of New Generations', by the Netherlands Organi- 
zation for Scientific Research (NWO/STW VIDI), and 
by the German Science Foundation (DFG) through its 
priority program 'Micro- and nanofluidics'. We acknowl- 
edge computing resources from the Jiilich Supercomput- 
ing Center and the Scientific Supercomputing Center 
Karlsruhe. 

Appendix A: Effective slip for a flow past strongly 
slipping anisotropic patterns 
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FIG. 6: The normal velocity gradients for the no-slip tex- 
tures with large amplitudes of slip length variation, bo/L = 1 
(dashed curve), bo/L — 5 (dash-dotted curve), bo/L = 20 
(dotted curve) as functions of stretched coordinates. 



large slope, of the order of L/r 3> 1 or larger, far from 
it. 



VI. CONCLUSION 

We have investigated shear flow past a super- 
hydrophobic surface with a cosine variation of the local 



In this Appendix we give some simple arguments show- 
ing that the flow becomes isotropic in case of large local 
slip, /3min = niin [/3 (y)] 3> 1. In this situation for arbi- 
trary local slip length /? (y) the leading-order bofi can be 
obtained directly from the local slip boundary condition, 
Eq.(|6]), without a need to solve the Stokes equations. In- 
deed, if ;9 > 1, Eq. ^ requires that the velocity at the 
surface is large, 

|Uslip + Ui^ {X,y,0) \ ^ /3min > 1, 

but the normal gradient is finite, \duiT {x, y, 0) /dz\ ^ 1. 
This can be fulfilled if only |usiip| ~ /3min 3> 1, |uit-| ^ 
1. One can then divide Eq. ([6]) by /3 to obtain 

slip duir I Q-\ \ 



z = 



u. 



/3 (y) dz 



Finally, by averaging of the last equation over the texture 
periods, and by using Eqs.(|4]) and Eq.([9]) we can easily 
derive an expression for a leading-order effective slip: 



L 



dy 

Hy) 



(Al) 



It can be verified that our result, Eq.(28), is consistent 
with (All. Note that the expression 



^eff 



^off 



1 - , 



(A2) 



8 



which has been obtained for stripes with 60, bi ^ i 
and is similar to the addition rule for resistors in parallel, 



also satisfies Eq.( Al ). 
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